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Abstract 



The Next to Minimal Supersymmetric Standard Model (NMSSM), proposed as 
a solution of the \i problem of the Minimal Supersymmetric Standard Model, has 
a discrete Z 3 symmetry which is spontaneously broken at the electroweak phase 
transition, resulting in a cosmological domain wall problem. In most cases this 
domain wall problem cannot be solved by explicit Z 3 breaking without introducing 
supergravity tadpole corrections which destabilize the weak scale hierarchy. Here we 
consider the possibility of solving the domain wall problem of the NMSSM via spon- 
taneous discrete symmetry breaking occuring during inflation. For the case where 
the discrete symmetry breaking field has renormalizible couplings to the NMSSM 
fields, we find that the couplings must be less than 10~ 5 if the reheating temperature 
is larger than 10 7 GeV, but can be up to 10 -3 for reheating temperatures of the order 
of the electroweak phase transition temperature. For the case of non-renormalizible 
couplings, we present a model which can solve the domain wall problem for large 
reheating temperatures without requiring any very small coupling constants. In this 
model the domain walls are eliminated by a pressure coming from their interaction 
with a coherently oscillating scalar field whose phase is fixed during inflation. This 
oscillating scalar field typically decays after the electroweak phase transition but 
before nucleosynthesis, leaving no additional Z 3 symmetry breaking in the zero- 
temperature theory. 



1. Introduction. 

The Minimal Supersymmetric extension of the Standard Model (MSSM) has 
become widely regarded as the most probable theory of physics beyond the Standard 
Model (SM), due to its ability to protect the weak scale from quadratically divergent 
radiative corections jl]]. However, the MSSM requires, in order to give an expectation 
value to both Higgs doublets, the introduction of a SUSY mass parameter for the 
Higgs doublets, u, which has no obvious reason to be of the order of niw, although 
some mechanisms have been suggested in the context of supergravity models 
(We will refer to this as the \i problem). One particularly simple and attractive 
solution of the u problem is to extend the MSSM to the Z 3 symmetric Next to 
Minimal Supersymmetric Standard Model (NMSSM) |3|, [|]. This adds a scalar 
superfield N to the MSSM whose expectation value provides an effective u term 
for the Higgs doublets. The u term is naturally of the order of the weak scale 
so long as the Z3 symmetry is imposed. However, the Z3 symmetry introduces a 
cosmological problem in the case where the reheating temperature after inflation is 
large compared with the weak scale, namely the formation of stable domain walls 
H, U due to sponanteous breaking of the Z 3 symmetry at the electroweak phase 
transition 0. 

In order to eliminate these domain walls, a Z3 symmetry breaking which has 
the same phase over the region corresponding to the observable Universe must be 
introduced into the effective weak scale theory. However, any explicit Z 3 symme- 
try breaking (even if due to Planck-scale suppressed non-renormalizible operators) 
will almost certainly destabilize the weak scale by allowing quadratically divergent 
gauge singlet tadpole diagrams which appear in the full supergravity theory |7|, §]. 
(However, for the special case where the Z3 is embedded in a gauged R-symmetry, 
such corrections can be avoided ||). 

The obvious alternative is to consider spontaneous Z 3 symmetry breaking due to 
a new gauge singlet scalar field, S, whose initial value is fixed during inflation. This 
ensures that the resulting discrete symmetry breaking phase is the same over the 
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whole observable Universe. In this paper we will discuss whether such a possibility 
can be implemented in a natural way, in the sense of not requiring extremely small 
coupling constants. We will consider two different approaches. One approach is 
based on coupling S directly to the NMSSM via renormalizible couplings. The 
second is based on having only non-renormalizible couplings of S to the NMSSM 
sector, suppressed by powers of the Planck scale or some intermediate mass scale. 

We first outline the cosmological scenario we have in mind. We will consider 
the simplest scenario, in which there is an initial period of inflation due to the 
energy density of an inflaton field, which subsequently oscillates coherently about the 
minimum of its potential until it decays, leaving the Universe radiation dominated 
at a reheating temperature Tr [|l(|. Tr should be low enough not to thermally 
regenerate gravitinos [|TT] , p| (Tr < 10 8 — 10 9 GeV for gravitino masses in the range 
lOOGeV to 500GeV (n|), implying that the value of the expansion rate of the 
Universe at Tr, H(Tr), should be less than lGeV. Between the end of inflation 
and the beginning of radiation domination the Universe is dominated by the energy 
density of the coherent inflaton oscillations. An important point to note is that 
the Universe does not "reheat" as such. During the inflaton oscillation dominated 
period the decay of the inflaton results in a themal background of particles with a 
temperature T r « k r {MpiHT R ^) 1 / 4 (where k r = (^ ^j ) ^ is approximately equal 
to 0.4 for temperatures larger than 100 GeV), the energy density of which eventually 
comes to dominate the energy density of the Universe at Tr ||10|| . The magnitude of 
the density perturbations observed by the Cosmic Background Explorer implies that 
H fa 10 14 GeV during inflation, assuming that the density perturbations are due to 
quantum fluctuations []14 |. Therefore it is possible that T r could exceed 10 9 GeV at 
the end of inflation, even if Tr is less than 10 9 GeV. In general T r is given by 

TrS36 . 3xl0 » fc ,.(_) ( w ^) GeV (1). 

Although it might seem that T r > 10 9 GeV would require a tighter bound on Tr 
in order to evade the thermal regeneration of gravitinos, it is straightforward to 
show that the rapid increase of if as a function of T during the inflaton oscillation 



dominated period, if oc T 4 , ensures that the gravitinos generated thermally during 
this period never exceed those generated around T R . Thus it is sufficient to impose 
Tr "Z 10 9 GeV. On the other hand, we will find that the higher temperatures which 
exist during the inflaton oscillation dominated period must be taken into account 
when considering constraints on the couplings following from the requirement that 
the Z% symmetry is not thermally restored before the NMSSM domain walls form. 

It has become clear that in many inflation models based on supergravity, the 
scalar fields will obtain order if 2 SUSY breaking mass squared terms as a result of 
the non-zero energy density which exists in the early Universe |10|, |i~5f . (A possible 
exception is the case of D-term inflation, depending on the details of the model 
HI). If such mass squared terms happen to be negative, then any scalar field with 



a flat scalar potential will have a large initial expectation value which is fixed during 
inflation. This is in contrast with the case of potentials having hidden sector SUSY 
breaking mass squared terms, which have a fixed value, m 2 , which is typically of the 
order of (100Ge\^) 2 . (By "hidden sector" mass terms we mean those coming from a 
supergravity hidden sector breaking mechanism In this case the fields cannot 
roll to the minimum of their potentials during inflation since H is much larger than 
m s and the fields are overdamped. Therefore they will have random initial values. 
The role of the if 2 mass terms in our discussion will be to fix the initial values of 
the fields at the end of inflation to be at a non-zero minimum of the scalar potential, 
with a single discrete symmetry breaking phase througout the observed region of 
the Universe. It is interesting to ask whether the domain wall problem can be solved 
with random initial values of the fields at the end of inflation, as in the case without 
the order if 2 mass terms, but we will not address this question here. 

There are then two possible scenarios for the cosmological evolution of the S 
expectation value, depending on the sign of m 2 . If this is negative, the 5* expectation 
value will remain large and constant once the hidden sector mass squared term 
comes to dominate the order if 2 term during the expansion of the Universe. On 
the other hand, if it is positive, the S scalar will begin to oscillate coherently once 
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m 2 s dominates the if 2 term. If these Z% symmetry breaking expectation values can 
survive down to temperatures less than that of the electroweak phase transition, T ew , 
(at which the Z 3 domain walls form), then we may be able to use the interaction of 
< S > with the NMSSM sector fields to eliminate the domain walls. 

However, we will see that it is not so easy to achieve a natural solution of the 
domain wall problem. On the one hand, we have to provide a large enough vac- 
uum energy density splitting between the different Z% domains to provide sufficient 
pressure to drive away the domain walls. On the other hand, for the case of the 
negative hidden sector mass squared term, we have to ensure that the Z 3 symmetry 
is not thermally restored while at the same time ensuring that the S expectation 
value does not give a large mass to the scalars in the NMSSM sector responsible for 
electroweak symmetry breaking. For the case of the coherently oscillating S scalar, 
we have to ensure that the S field is not thermalized before the electroweak phase 
transition and that it decays after the electroweak phase transition but not so late 



that it disrupts the predictions of standard big-bang nucleosynthesis |13 |. 
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2. Renormalizible S couplings to the NMSSM sector. 

(a) m 2 < 0. 



We define the transformation of a field with charge Q(<f>) under a discrete 
symmetry by (f> — > e l27r< ^</> (— \ < Q < |). The simplest extension of the NMSSM is 
then obtained by considering the Z 3 charge of S to be equal to that of N, Q(S) = 
1/3. The superpotential of the model is given by W = Wnmssm + Ws, where 

k 

Wnmssm = Wyukawa + XnNHiH 2 — —N (2) 

is the standard Z 3 symmetric NMSSM superpotential |], |3|, £|, 0] and 

W s = XshSH^ + ^-S 3 + X 2 S 2 N + X 3 SN 2 (3) 

gives the couplings of S to the NMSSM fields. In addition we will have soft SUSY 
breaking terms of the form [ITS 



(m2±# 2 )|0,| 2 + (A Q ^ Q + /i.c.) (4), 

where <pi are the scalar fields and W a are the trilinear superpotential terms. (A a ~ 
m s ±H in the presence of non-zero H, with the sign of the H correction being model 
dependent [0). With these superpotential and soft SUSY breaking terms, the S 
scalar potential is assumed to be given by 

V(S) = (-m 2 s - H 2 )\S\ 2 + (A X ±S 3 + h.c.) + (A 2 + A 2 )|S| 4 (5). 

The expectation value of S is then given by 

, q ^„ (ra 2 + H 2 ) 1/2 

<b> ~ {xi + xiy 2 {b} - 

(We have neglected the A Xl term, which can alter < S > by at most factor of the 
order of 1). 

In this we have ignored the N expectation value. In fact, a non zero < N > will 
generally exist at the minimum together with < S >. This could result in a large 
effective fi term and affect the natural generation of the weak scale. We will discuss 
iV expectation value and its consequences later in this section. 



We first consider the conditions under which there is no thermal restoration of 
the Z 3 symmetry. This requires that the S 2 term in the finite temperature effective 



potential |12], [17j remains negative for all values of T for which the S condensate 
scalars are in thermal equilibrium. 

The iV particles and Higgs particles Hi will contribute to the S finite-temperature 
effective potential if they are lighter than T. In fact, they will be lighter than T for 
T large compared with T ew , since any N or Hi mass term due to < S > cannot be 
very large compared with the weak scale without disrupting electroweak symmetry 
breaking ||. The S particles orthogonal to < S >, which we denote by S' , will 
also gain a mass from < S > and can contribute to the finite temperature effective 
potential. The condition for these particles to be light compared with T for H < m s 
is that 

Al < T (7), 



which will generally be satisfied for T > T ew . (We need not consider the constraints 
for H > m s , since in this case the soft SUSY breaking mass squared in the scalar 
potential becomes of the order of H 2 , which will increase with T much faster then 
the thermal T 2 correction, thus ensuring that the Z 3 symmetry remains broken at 
H > m s if it is broken for all H ^ m s ). The effective SUSY iV mass from the A3 
coupling is then given by 7tin{S) = 2X3S. Similar contributions will come from Ai 
and A 2 . 

We must then ensure that the temperature dependent contribution to the S 2 
term in the finite temperature effective potential due to mjv, m%(S)T 2 /8, is small 
compared with m 2 S 2 for all T ^ T eq , where T eq is the temperature at which the S 
scalars come into thermal equilibrium. The dominant process bringing the S scalars 
into thermal equilibrium with the light particles in the thermal background will 
be inverse decays, with a rate r m „ fa ^^-T (where the factor is not too small 
compared with 1; for the case of stop quarks this has been estimated to be typically 
of the order of 0.1 fl8|l ). If T eq < T R , then the thermal equilibrium condition T inv ^ H 
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occurs during radiation domination and T eq is given by 

where the expansion rate during radiation domination is given by H — fcr j^ T , with 



kx(T) « 16 for T larger than around lOOGeV [0. This is consistent with T eq < Tr 

If this is not satisfied then T eq > Tr and, with H(T) given by H(T) = T A (kjT^,Mpi) l , 
the thermal equilibrium temperature becomes 

^-(^yv^M^ (io). 

The Z 3 symmetry remains broken at T e? if Aj ^ m s /T eq . For T e? < Tr this requires 
that 



where we have used ~ 16. This is the condition for avoiding Z 3 symmetry 
restoration if Tr > 10 7 GeV A . If Tr is less than 10 7 GeV, the upper limit on Aj occurs 
when T e(J > Tr and the condition that the Z3 symmetry remains broken at T eq 
becomes 

A-2.1X10- K /^— ) (^) (12). 

For T fi = 10 2 Ge1/ (10 5 Gey) this requires that A, < 10~ 3 (10~ 4 ). Thus we see that, 
for the case of m? s < and renormalizible S couplings to the NMSSM sector, the 
couplings of S to N and to the Higgs fields must be rather small if the reheating 
temperature is large compared with T ew , in order not to thermally restore the Z 3 
symmetry. In particular, if the reheating temperature is close to the gravitino upper 
limit, then these couplings must be less than around 10 -5 . 

We have not yet discussed the elimination of the domain walls. We will discuss 
this in more detail later, but for now we merely note that it is very easy to introduce 
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large < S > dependent contributions into the N scalar potential. For example, from 
the superpotential couplings k and A2, we obtain a term in the scalar potential 

-2\ 2 k < & > 2 N 2 + h.c. ~ (^ Xf + xl j km2 s( N2 + h - c -) ( 13 )- 

However, if we assume that k is not very small compared with 1, such that 2k > X 2 , 
then A2 must be smaller than \\/k in order to avoid introducing a mass term for the 
N scalars much larger than \Q 2 GeV . In the NMSSM broken by radiative corrections, 
the phenomenologically favoured region of parameter space has < N > > ITeV. In 
this case the N particles essentially decouple and one has an effective MSSM with 
= Atv < N > and a contribution to the soft SUSY breaking "j3-term" given by 
v = k < N > [|J. Both of these cannot be very much larger than 10 2 GeV without 
disrupting the natural generation of the weak scale. (It is possible to have very 
large < iV >, so long as A a? and k are correspondingly small. The condition for a 
weak scale /1 term is then that A k ~ 1 ) . The effect of a large mass term from 
equation (13) would be to introduce an effective negative mass squared, m 2 N for 
the N scalars, with the N expectation being given by < iV >~ \rriN eff\/k, which 
would result in v ~ |mjv e //| @j- Requiring that v be small compared with a mass 
scale m* not much larger than 10 2 GeU and that the no-thermalization condition 
(12) is satisfied then implies that 

d \ k )\ T R ) VlOOGe^y KimGeV) y ! 
(For Tr ^ itfGeV the upper bound from (14) with Tr ~ itfGeV applies). Thus for 
large reheating temperatures, Tr > 10 7 GeV r , we would require that A2 "Z 10 _7 (0.01//c) 
in order to avoid a mass term for the N scalars much larger than lQ 2 GeV . For smaller 
Tr this upper bound becomes weaker, allowing A2 to be as large as 10 _3 (0.01//c) 
for Tr 10 2 GeV . Therefore, unless k ^ 10~ 4 (and, since X^/k ~ 1 is necessary 
[f|], Aat ^ 10~ 4 ), for large Tr there will be a tighter upper bound on some of the 
couplings than that coming from Z3 symmetry restoration. 

We next consider the N expectation value. When < S > is introduced into the 
N — S scalar potential, terms linear, quadratic and cubic in N arise, which result 
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in a non-zero < N >. In general it is difficult to discuss the minimization of the 
potential analytically, but for the natural case where the couplings A, (i=l,2,3) all 
have the same magnitude and \k\ is large compared with |Aj|, and where < N > 
does not alter the value of < S > significantly, we find that < N > is essentially 
determined by the linear and quartic terms in the N scalar potential. In this case 

< iV > is given by 

/\\ 2 / 3 

< N>&lj\ <S> (15). 

(This will not affect the value of < S > so long as < iV > is small compared with 

< S >). The main effect of < N > will be to introduce an effective /i term for the 
Higgs doublets via the A a? superpotential coupling, 

(r) ms (16) - 

Since A^ and k must be of the same order of magnitude to naturally generate a weak 
scale \i term, we see that A, cannot be too small compared with k, thus favouring 
large A, and/or small k. We note that < S > also introduces a n term via the Xsh 
superpotential coupling. This will be acceptable so long as Xsh is not larger than 
(A? + Ai)V2. 

We should check that spontaneous breaking of the discrete symmetry does not 
introduce dangerous quadratic divergences. A quadratically divergent tadpole arises 
at two-loops in N = 1 supergravity with hidden sector SUSY breaking and no direct 
coupling of the hidden and observable sectors in the Kahler potential ||. With a 
Kahler potential given by 

K = N^N + HlH l + o[^ (17), 

corresponding to minimal scalar kinetic terms, and superpotential terms of the form 



W = ... + XNH 1 H 2 - -^N 2 ^^ + O ( -L 

Mpi \M Pl/ 

a divergent iV tadpole is generated 



^/AW^±^ (19). 
(16tt 2 ) 2 7 \ M Pl J 1 ; 
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On introducing the hidden sector SUSY breaking F-term for K, Ff c 9 2 9 with Fk 
( "v^' ) ' ^ nere * s a contribution to the A" scalar potential of the form || 



2AA a A 2 (m s M Pl \ 2 



(N + N ] ) (20). 



(16tt 2 ) 2 M^ ^ y^F , 

Taking the cut-off A to be Mpi, we find that the additional term in the A" scalar 
potential is of the form M 3 (N + Aft), where M « 7x 10 5 X 1/3 X 1 J 3 GeV . If M is much 
larger than 10 2 GeV A , this term will make impossible the generation of a naturally 
weak scale /x term. Throughout this paper we will consider the natural value of 
non-renormalizible couplings such as X a to be of the order of 1, with the strength 
of the couplings set by the large mass scale which we will take to be the Planck 
scale. In this case we would require that A < 10 -12 in order to have M ^ lOOGeV^. 
In the case where discrete symmetry breaking is spontaneous and non- vanishing at 
zero temperature, we can still have quadratically divergent tadpoles. For example, 
for S with charge Q(S) = 1/3, one can have a non-renormalizible operator of the 
form of the X a term but with an overall factor of (jjj^f) ■ F° r S with charge 
Q(S) = 1/6 this factor becomes (-f^ ) • As a result, the quadratically divergence 
is not a problem for all A < 1 if < S > < 10 13 CreV for the Q(S) = 1/3 case and 
< S > < 10 16 GeV for the Q(S) = 1/6 case. From equation (6) we see that, for 
m s « 10 2 Ge1/, < S > < 10 13 GeV requires that (A 2 + A 2 ) 1 / 2 > 10 -11 . This can be 
easily satisfied. 
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K<J^\(^) (2D, 



(b) m| > . 

In this case we will have a coherently oscillating S scalar. We have to ensure 
that the Z 3 -breaking condensate does not thermalize or decay before the electroweak 
phase transition. 

The requirement that the condensate does not thermalize via inverse decays 
imposes the constraint, 

1/2 j, N 1/2 

K d ) \~Mpi. 

where we have assumed that the Universe is radiation dominated. Thus, with T w 
T ew , this requires that Aj < 10 -8 . (Below T ew the S' , Hi and N particles gain 
masses from electroweak symmetry breaking). This constraint should be applied if 
the masses of the S' , Hi and N particles due to < S > are small compared with T. 
During the inflaton oscillation dominated era, the energy density of the S oscillations 
scales as that in the inflaton oscillations, whilst after reheating the energy density 
of S oscillations scales as a(T)~ 3 , where a(T) is the scale factor. Thus the energy 
density during radiation domination at T < Tr is given by 

_ a p (T)T 3 T R 
Ps{T) - mlMl, Po (22) ' 

where p a is the initial energy density of the S condensate and ot p (T) = 4^ kr (Tr) 2 . 
The the corresponding amplitude of the S oscillations at temperature T is S(T) w 

1 /2 

(Ps/Po) S , where S a is the initial amplitude of the oscillations. The condition 



that the mass of the particles due to their interaction with the condensate [jlS 
XiS(T), is less than T then becomes 

Xt - [-^rfT- R ) (23) - 

This will be easily satisfied for T ew < T < Tr. Therefore (21) applies. We should 
also check that the energy density is not dominated by the coherent S oscillations. 
This is satisfied if + A| > 10~ 28 m 2 s /T, which is typically easily satisfied. 

Under the assumption that the S condensate particles are heavy enough to decay 

A 2 

directly to a two-body final state via the coupling Aj, with a rate Td w ^m s , the 
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requirement that the condensate does not decay before temperature T imposes the 
constraint 



This will give an even tighter constraint on Aj than the no-thermalization constraint 
(21) if the coherent oscillations have to survive down to a temperature much lower 
than the electroweak phase transition in order to eliminate the domain walls, as is 
likely to be the case in the NMSSM |§. 

We next consider the conditions under which the S condensate can eliminate the 
NMSSM domain walls. This requires that the pressure due to the energy difference 
AV between the domains is large enough to cause the higher energy domains to 
collapse on the time scale H{T)~ l || ||, |20fl . The condition on AV depends on 
whether the domain walls are relativistic or not, i.e. whether there are significant 
frictional forces due to the interaction of the domain walls with the thermal back- 
ground particles. If there are strong frictional forces, then the walls will expand 
until the force due to surface tension is balanced by the friction, which will occur 
at a radius much less than the horizon H~ l (T). The pressure due to AV must 
then overcome this large frictional force. If there are effectively no frictional forces, 
meaning that the surface tension can overcome the frictional forces for relativistic 
walls of radius equal to the horizon, then the condition on AV will be typically much 
weaker. For the case of an unbroken Z 3 symmetry, the conventional view would be 
that there are no (or at least highly suppressed) frictional forces, since the particles 
crossing the domain walls would have equal masses on either side of the wall and so 
would transfer no net momentum to the wall on passing through |2(], |2l|. However, 
it has been suggested that, on taking into account the change in the particle masses 
in the vicinity of the wall, there will be a non-zero reflection coefficient, leading 
to a large frictional force until the heavier particles become Boltzmann suppressed 
at temperatures less than around that of the quark-hadron phase transition, after 
which the domain walls become relativistic on the scale of the horizon 0. In more 

5 

detail, the friction force per unit area on the wall is estimated to be f ~ S^i^rvT , 




(24). 



13 



where m is the mass of the heaviest thermal background particle with m < T, g is 
the number of degrees of freedom of the particle and v is the velocity of the wall [?]]. 
This force rapidly decreases as T decreases and the heavier thermal particles be- 
come Boltzmann suppressed. The condition on AV to eliminate the domain walls, 
AV > /, will be most easily satisfied once the frictional force is weak enough to al- 
low the domain walls to remain relativistic up to scales of the order of the horizon. 
Thereafter, the condition will become more difficult to satisfy as T decreases, since 
the pressure difference between the domains rapidly decreases as < S > decreases 
with T. 

Suppose the domain walls become relativistic on horizon scales at T re [. Soon 
after T re \ , the uncollapsed domain walls will have a radius of the order of the horizon 
H(T re i)~ l . This will continue to be true until the pressure due to AV dominates 
the force due to the surface tension a. This in general occurs once the radius of a 
domain is larger than a critical radius, r c , given by 

(25). 
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The smallest |AV| will therefore correspond to the largest possible r c , the horizon 
radius. Although the Z 3 domain wall will consist of the Standard Model Higgs fields 
as well as the N field, for simplicity we will consider the wall to be made only of 
the N field. In this case the domain wall potential will be due to the soft SUSY 
breaking N 3 A-term and the surface tension will be o ps m s < N > 2 H [2CJ . 

From the superpotential couplings A2 and k we obtain an energy density splitting 
between the different Z 3 vacuum states 

AV ps 2X 2 kS 2 N 2 (26). 

This depends on the value of S 2 and so will have a non-zero value on averaging 
over the S oscillations. The condition that this can overcome the surface tension 
for horizon-sized relativistic domains at T re i is then that 



^1 + ^2 < a p(Trel) TrT, 



2kX 2 " k T (T re i) m s Mpi 
14 



'" ; (27). 



The left-hand side will be smallest for Ai < A 2 , so we require that 

Ai < A2 < 1Q- I3 a,(r ref ) / T R \ / T rel \ /lOOGe^X 

fc ~ fc ~ A; T (T re/ ) U0 9 GeVV U.lGeW V m s / 1 J ' 

Therefore we see that the requirement that AV is large enough to eliminate the 
domain walls imposes a very tight constraint on at least some of the couplings. 
A similar S 2 dependent contribution to AV, AV ~ 2AiA 3 5' 2 iV 2 , arises from the 
superpotential couplings Ai and A 3 . For this to eliminate the domain walls we 
would require that 



A 



2 < Ai < lQ- 13 a p (T rel ) / T R \ f T rel \ f imGeV 



A 3 ~ A 3 ~ k T {T rel ) \WGeVJ \0.1GeVj 

which would demand even smaller values of Ai and A2 once the no-thermalization 
and decay constraints, (21) and (24), are imposed. 

Thus we can conclude that, for the case of renormalizible couplings of the Z 3 
symmetry breaking scalar S to the NMSSM, very small couplings (less than 10~ 8 
and in some cases less than 1CT 13 ) are required in order to eliminate the NMSSM 
domain walls if the hidden sector mass squared of S is positive. If it is negative, small 
couplings < 10~ 5 (and even smaller if the NMSSM couplings k and A^ are larger 
than 10~ 4 ) are again required if the reheating temperature is larger than 10 7 GeV, 
but, for smaller values of Tr, the constraints from avoiding thermal restoration of 
the Z 3 symmetry become somewhat weaker, allowing couplings as large as 1CT 3 for 
Tr ~ T ew . 
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3. Non-Renormalizible S couplings to the NMSSM Sector. 

We have seen that the case of renormalizible S couplings to the NMSSM sector 
requires rather small couplings when the reheating temperature is large (less than 
1CT 5 for the case where the reheating temperature is larger than 10 7 GeV, with 
an even smaller upper bound if the NMSSM N couplings are not less than about 
10 -4 ). In particular, for the case of a positive hidden sector mass squared for the 
S scalar, some couplings must be less than about 1CT 13 , regardless of the reheating 
temperature. One way to overcome the need for such small couplings might be 
to simply eliminate all the renormalizible couplings of the S scalar to the NMSSM 
fields, allowing only non-renormalizible couplings suppressed by powers of some large 
mass scale (which we will take to be the Planck scale). This is naturally achieved 
by extending the Z 3 symmetry of the NMSSM to a discrete symmetry, which we 
will call Za, under which the charge of S is such that no renormalizible couplings 
of S to the NMSSM fields are possible, 
(a) m 2 < 0. 

In this case there will be a large and constant expectation value for the S scalar 
once H < m s . This large expectation value will introduce what we believe to be a 
generic problem in this case, namely a very large SUSY mass for N, much larger 
than the weak scale. To see this, consider the simplest example, corresponding to S 
having a Za charge 1/6 whilst N, as usual, has charge 1/3. Then S can only have 
leading order non-renormalizible superpotential couplings of the form 



M P , M P , 12 6M? 



Lpi IVlpi U-iWpj 

The expectation value following from the S scalar potential, 

V(S) = -m 2 s S 2 + (31), 

gives, with m s pz 10 2 GeV and As « 1, < 5 >^ 5 x 10 14 GeV. This results in an 
SUSY mass for N of the order of 10 10 GeU, making a solution of the \x problem 
impossible in this case. We find that decreasing the Z A charge of S or introducing 
an intermediate mass scale (in the manner of the model to be described in the 
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following section) only increases the N mass. Based on this, we conclude that the 
case of S scalars with non-renormalizible couplings and a negative hidden sector 
mass squared is unlikely to provide a natural solution of the NMSSM domain wall 
problem. 
(b) m\ > . 

So far we have found no model which can solve the NMSSM domain wall prob- 
lem for reheating temperatures much larger than the electroweak phase transition 
temperature without requiring small couplings. We now present a model, based on 
a coherently oscillating scalar, which can eliminate the domain wall problem of the 
NMSSM for values of Tr as large as the gravitino upper bound and which does not 
require any small couplings. 

The model we consider requires, in addition to the NMSSM singlet N and the 
Za breaking field S, two additional singlet fields, X and B. X will acquire an 
intermediate mass as a result of an expectation value for B. The model is also 
assumed to have a more complicated discrete symmetry, Za X Zb, where Za is the 
extension of the original NMSSM Z 3 whilst Zb is a discrete symmetry under which 
all the NMSSM fields are singlets. Zb is introduced in order to control the non- 
renormalizible terms. The discrete symmetry charges of the fields are as defined in 
Table 1(a). 

Table 1(a). Field charges under Za x Z b - 





S 


X 


B 


N 


Z a 


1 

6 


1 

2 





i 

3 


Z b 


1 
3 


1 

3 


1 

6 






It may well be significant that the Za x Zb charges of the remaining NMSSM 
fields can always be chosen to eliminate the dangerous renormalizible B and L 
violating operators u c d c d c , d c QL and e c LL 0. For example, this is true for the 
charges of Table 1(b). Therefore the absence of these renormalizible interactions 
could be interpreted as a sign of the existence of a non-trivial discrete symmetry 
such as Za x Zb, rather than of a simple R-parity or conservation of B and L. 
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Table 1(b). Example of field charges eliminating dangerous renormal- 
izible operators. 





e c 


L 


d c 


H d 


Q 


u c 


H u 


Z A 
Zb 


1 

6 

1 

3 


i 

2 

1 

3 


1 

3 




i 

3 

















The superpotential of the model is then given by 

W = W N mssm + W x + W NR (32), 
where Wx is the allowed renormalizible coupling of X, 

W x = X a XSN (33) 

and Wnr gives the leading order non-renormalizible superpotential couplings in- 
volving the additional gauge singlet scalars, 

W NR = 4§rS 6 + ^B 2 X 2 + -4§tB 6 (34). 
QM Pl Mpi 6M Pl v ; 

The model will require an intermediate mass for the X scalar. This will be generated 
by an expectation value for B. The obvious possibility is to assume that the hidden 
sector mass squared of the B scalar is negative, which would result in an expectation 
value for B given by, 

^5xl0 14 GeV (35). 

However, we should perhaps note that in this case there could be problem due to 
the late decay of oscillations of the weakly coupled B scalar around the minimum 
of its potential. The B scalar has a mass of the order of m s , whilst the X particles 
to which it couples gain a mass of the order of lO^GeV as a result of < B >. 
Thus the B decay rate will be highly suppressed in this case. If the B field is 
not very close to the minimum of its potential at the end of inflation, its coherent 
oscillations about the minimum of its potential will decay long after nucleosynthesis 
and could dominate the energy density of the Universe when they decay. This is 
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< B > 




all dependent upon the value of B at the end of inflation and the dynamics of its 
subsequent evolution as the Universe expands and its effective (negative) 0(H 2 ) 
SUSY breaking mass term decreases. Since in our model the B field only serves to 
generate an intermediate mass for X and, in a sense, may be regarded as a toy field 
representing the dynamics of a more realistic model, we will not concern ourselves 
with the details of this issue here. 

In general we will consider Mx ^ 10 10 GeU. One reason for this choice, as we will 
show later, is that it ensures that the X particles are Boltzmann suppressed for all 
temperatures, thus eliminating the possibility of thermalization of the S condensate 
by light thermal X particles. Integrating out the massive X field, we obtain the 
following effective interaction in the superpotential 

^-S 2 N 2 (36). 
M x 

This is valid so long as A a < S > < M x , such that the mixing between the X and 
N fields via the coupling (33) may be neglected. This interaction is fundamentally 
important in what follows, as it will allow both a sufficiently rapid decay of the S 
oscillations so as to avoid problems with nucleosynthesis and at the same time allow 
a sufficiently strong interaction of the S oscillations with the N scalars to drive away 
the Z 3 domain walls. Assuming that the S scalar mass is large enough, the S scalars 
will decay to iV particles and S fermions (which will be highly decoupled and will 
typically have a mass not much larger than leV, coming from the interaction (36) 
with < iV >~ lOOGeU — ITeV). So long as their energy density is sufficiently small 
compared with than that of the radiation when they are produced, the S fermions 
will not give rise to any cosmological problems. (This may be difficult to achieve, 
however, if the Universe is dominated by the S scalars when they decay, since one 
would expect approximately 1/3 of the energy density from S decays to be in the 
form of S fermions. Alternatively, for the charge assignment of Table 1(b), one can 
obtain an operator of the same form as equation (36) but with S 2 N 2 — > SNH U L, 
which can allow the S scalars to safely decay to NMSSM fields). 
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The S scalar potential is given by 



V(S) « (m 2 s - H 2 )S 2 + i^S 5 J (37). 

The initial value of the S field when the oscillations begin at H « m s , S^, is therefore 
given by 

We first consider the conditions under which the S energy density decays at a 
temperature less than that of electroweak phase transition, T ew , without disturbing 
the predictions of nucleosynthesis. This requires either that the ratio of the energy 
density in the S scalars to that in radiation when the S scalars decay, r D , is less 
than 10~ 6 , in order not to photodissociate the helium abundance O, E|, or that 



2 



the S scalars decay at a temperature greater than around lOMeV [p3| . The former 
condition requires that 

T R <^MkT d (39). 

The S scalars in the condensate decay via the interaction of equation (36) with a 
rate given approximately by 

a d X 4 a m 3 s 



Ml 



(40), 



where ~ 1927r3 ~ 10~ 4 and we are assuming that the S scalars are heavy compared 
with the N particles. This gives for the decay temperature, assuming that decay 
occurs during radiation domination, 

Therefore, in the case where decay occurs after nucleosynthesis at Td ^ lMeV, we 
require from equation (39) that the reheating temperature satisfies, 

T R <r D m sl \l(^l) (42), 



x 



where 

7 ^k T (T d yn 
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Thus with r D « 10" 6 , 7 « 10~ 4 , m 5 « 10 2 GeU and M x > 10 10 GeU, this requires 
that Tr < a 1 J 2 Xll0 4 GeV . This imposes a severe restriction on T R for A a < 0.1. On 
the other hand, from equation (41), we see that it is quite likely that T d > lOMeV. 
In general, T d is given by 

T - K 10 lio^) (oij liooG^J (^rj (43) - 

Therefore, in this case we expect that, for M x ~ 10 10 GeU (which will be seen to 
be the preferred order of magnitude for M x in this model), the S oscillations will 
typically decay at a temperature between lOMeV and lGeV for A a in the range 0.1 
to 1. The amplitude of the S oscillations at T < Tr (assuming that the Universe is 
radiation dominated) is given by 

16a p (T)V2/ ms \ 3/4 , T , 3/2 T 1/2 

<5> ^^^(too^v7) Uj (low) Gey (44) - 

The corresponding energy density in the coherent 5 oscillations is ps ~ m 2 s S 2 . This 
will come to dominate the radiation energy density once T ^ Td om , where 

r - (lo^f ( w) ^ < 45 '^ 

Thus if the domain walls become relativistic below T dom , the elimination of the 
domain walls by the condensate will occur whilst the Universe is dominated by the 
energy density of the oscillating S field and possibly whilst the background radiation 
is dominated by radiation coming from the decays of the condensate scalars. In 
the case where the S scalars decay during matter domination, it would seem that 
equation (41), derived assuming radiation domination, is not correct. However, since 

1 /2 

the decay condition is given by T d m H oc pg and since the radiation energy of the 
Universe when all the S scalars decay will simply be given by ps at this time, the 
decay temperature will be the same as that given in equation (41). 

We must also ensure that the S condensate is not thermalized by scattering 
processes involving the thermal background particles. We first show that the X 
particles are always heavier than T and so Boltzmann suppressed if the coupling A a 
is not too small compared with 1. The mass of the X particles from the X — N 
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mass matrix is given by M x ~ Max(Mx, X a < S >). Therefore M x Z 10 10 GeV if 
M x t lO w GeV. Since, from equation (1), T r at H xs m s is given by 



TJH « m. 



we see that Mx ^ T for all if less than m s . At if ~ m s the X mass will be given 
by M x ~ A a < 5 > so long as A a > 2 x 10~ 5 Ai/ 4 ( 10 i!fg e y ) • Assuming that this 
is satisfied, we see that since < S >oc H oc T 4 for if > m s , Mx will be larger 
than T for all temperatures after the end of inflation. Thus the X particles will not 
thermalize the condensate for A a > 1CT 5 . 

We next consider whether the N particles can thermalize the condensate. The 
N particle mass from the X — N mass martix is given by M N Min(X a < S > 
, Xa ^f^ )■ In the following, for simplicity, we will concentrate on the case of large 
Tr, not much smaller than the gravitino upper bound, which is of most interest 
to us here. Smaller values of Tr can be analysed in a similar way. Let T/v be the 
temperature below which Mjy becomes less than T. Let T c be the temperature 
below which \ a < S > < M X - Since for T > T c the N and X particles have the 
same mass and we have shown than Mx ^ T for all T, Mjy can only become less 
than T at temperatures less than or equal to T c . Assuming that T c < Tr, T c is given 
by 

Tc ~ A 2 / 3 UoGeW U^GeW ^ T R ) ^ (47) ' 

Therefore with A a not too small, say in the range 0.01 to 1, this will be less than 
Tr for Mx ~ 10 10 (jeV and for reheating temperatures not much smaller than the 
gravitino upper bound. M N can only be less than T if T N < T c . In this case T N is 
given by 

1.2xl0 5 Ag /4 / m s \ 3 / 4 /M x \ 1/2 



This is indeed less than T r if 
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This will be easily satisfied for values of Tr not much smaller than the gravitino 
upper bound, M x ~ 10 10 GeU and A a in the range 0.01 to 1. Therefore in this 
case we must apply the no thermalization condition during radiation domination 
at T/v- Thermalization will primarily occur via the effective interaction of equation 
(36). We treat the scattering of the thermal N particles from the zero momentum S 
condensate scalars as a simple scattering process involving S and iV particles. The 
scattering rate for SipN ^ SipN is then given by 

Requiring that the no thermalization condition, T sc < H, be satisfied at T/v, equation 
(48), then gives for the no thermalization constraint on Tr, 

TR > 6xl ^ xl (^Y' 2 (^L) 3 GeV (51) . 

H ~ °U00GeW \ M x J 

\ a < 0.2 and M x ~ 10 10 GeV will easily allow a range of values for T R below the 
gravitino upper bound to be consistent with no thermalization of the S condensate. 
Thus with, for example, A a « 0.1, M x ss \Qi w GeV and T R w 10 7 - g GeV, the 
S condensate will evade thermalization and will safely decay at a temperature of 
around lOMeV. 

We next consider the conditions under which the S condensate can eliminate the 
NMSSM domain walls. We first consider the case where the Universe is radiation 
dominated. The terms in the SUSY scalar potential responsible for the energy 
density splitting are given by 

^(S 2 NN^ 2 + h.c.) (52). 
M x 

This is explicitly dependent on the phase of the S field and has a non-zero average 
over time. It results in a splitting of energy density of the different N vacuum states, 
AV, given by 

AV w ^ < S 2 >< N > 3 (53). 
M x 

Thus the condition for AV to be large enough to eliminate relativistic domain walls 
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at a temperature T during radiation domination, AV(T) > aH(T), becomes 

Since < S 2 > is proportional to T 3 , we see that this is most easily satisfied as soon 

as the walls become relativistic at T re i. Thus in order to eliminate the domain walls 

during radiation domination we require that 

> k(T rel ) 1 m s m 2 s M x M Pl 
R ~ a p (T rel )k\lT rel <N> SI [ j ' 

For example, if we consider the domain wall to become relativistic shortly after 

the quark-hadron phase transition [7], say at T re i « O.lGeV, and to be too strongly 

damped for the pressure due to AV to be able to collapse the domain walls at higher 

temperatures, then we obtain 

T > s v m3^ /2 ( ™. \ ( m s \ 3/2 f Q.lGeV \ ( M x \ 



'-rel 

We have included a correction factor, f c , for the case where the domain walls become 
relativistic during S condensate matter domination rather than during radiation 
domination. We have shown that the Universe becomes matter dominated at a 
temperature Trf om , equation (45), which is in the tens of GeVs for Tr close to the 
gravitino upper bound. In this case, for T re \ < Td om , we have to modify the domain 
wall constraint on Tr. Once T ^ Td om , entropy will be effectively conserved and the 
radiation energy density will be mostly due to the primordial radiaiton. However, 
there will also be additional radiation due to the decay of the S condensate, with a 
temperature given by equation (1) with the replacement Tr — > This radiation 
will dominate the primordial radiation once T < T*, where 
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d) \ rpi/Brpl/5 



T^TjZ (57). 



Once T < entropy is no longer effectively conserved and H oc T 4 drops much more 
rapidly with T. The result of all this on the domain wall bound is to introduce a 
correction factor f c , where, for T re i > T, 

/ C =(^) V2 (58) 
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and for T re i < T„ 



j(T rel ) 



'rel 



The weakest lower bound on Tr is obtained when T rd > T*. T rei > T* is satisfied if 

rr, < ,1,8 f 1 \ V V ^re/ / WOGeV ^ f 10 9 GeV V f \ ^ r , , 

T ^ 40A " fed (tu^7) (— ) HH MeV (60) - 

The lower bound on Tr from requiring the elimination of the domain walls then 
becomes 

For example, this allows a range of reheating temperatures from Tr 5 x NftGeV 
up to the gravitino bound to be compatible with elimination of the domain walls at 
T rel « O.lGeV when A a is in the range 0.01 to 1, k « 0.1 and M x ~ 10 10 Ge\/. 

Thus we see that, with A a « A; « 0.1, M x « 10 10 GeV/ and T^ « 10 7 - 9 GeV/, 
the 5* condensate can evade thermalization, can safely decay before nucleosynthesis 
(at a temperature of around lOMeV) and can eliminate the NMSSM domain walls, 
without requiring any small renormalizible couplings. 

We should note that the preferred order of magnitude for M x in this model is 
around 10 10 GeV. Values smaller than this make it difficult to satisfy the no ther- 
malization constraint, equation (51), whilst values larger than this make it difficult 
to eliminate the domain walls, equation (61), and to ensure that the condensate 
decays before nucleosynthesis, equation (43). 

So far we have not explained why it was necessary to introduce the second 
discrete symmetry Zb- In the absence of this discrete symmetry, we could have had 
an additional term in the non-renormalizible superpotential, 

^S 3 X (62), 

Mpi 

where A& ~ 1. On integrating out the X field, this would give an effective operator 
in the superpotential, 

-.S'.V (63), 



M X Mpi 
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which lifts the flat direction much earlier than the S* 6 superpotential term, such that 

^^MMf* (64) 

As a result, the lower bound on Tr from the requirement that the domain walls can 
be eliminated during radiation domination would become 

Tr > 3 x 10 12 ^- ( ms V /3 ( 0AGeV ) ( M * ) 1/3 CcV (65) 

Even with k ~ A a ~ 1 this would require that Mx ^ 10 2 GeV in order to allow 
reheating temperatures below the gravitino bound, making it impossible to avoid 
thermalizing the condensate. 

Thus, with non-renormalizible couplings of S to the NMSSM fields, the Z 3 do- 
main wall problem of the NMSSM can be solved even if the renormalizible couplings 
are large and the reheating temperature is not small compared with the gravitino up- 
per bound. The introduction of an intermediate mass scale of the order of 10 10 GeV 
is essential for the solution of the domain wall problem in this case, since it provides 
non-renormalizible operators which are strong enough to enable the oscillating scalar 
field to decay fast enough to avoid problems with nucleosynthesis and to allow it to 
provide a sufficiently large pressure to eliminate the NMSSM domain walls whilst 
still being weak enough to prevent thermalization of the discrete symmetry breaking 
scalar condensate before the NMSSM domain walls have formed. 
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4. Conclusions. 

We have considered the possibility of solving the NMSSM Z 3 domain wall prob- 
lem by spontaneous discrete symmetry breaking occuring during inflation. For the 
case of renormalizible couplings of the discrete symmetry breaking scalar S to the 
NMSSM sector we find that, if the S scalar has a negative hidden sector mass 
squared term, then it is possible to solve the NMSSM domain wall problem without 
requiring extremely small couplings so long as the reheating temperature is not very 
large compared with that of the electroweak phase transition. Reheating temper- 
atures larger than 10 7 GeV would require the couplings to be less than 1CT 5 . A 
solution with a positive hidden sector mass squared for S always requires very small 
couplings and so seems generally unnatural. 

For the case of non-renormalizible, Planck-scale suppressed operators, a negative 
hidden sector mass squared solution appears generally ruled out by the generation of 
a large SUSY mass for the N field of the NMSSM. For the case of a positive hidden 
sector mass squared, we find that it is possible to solve the NMSSM domain wall 
problem via a coherently oscillating S scalar without any small couplings, even if 
the reheating temperature after inflation is large. This requires the introduction of 
an intermediate mass scale of the order of 10 10 GeV, in order to allow the S scalars 
to decay without cosmological problems, and an additional discrete symmetry in 
order to control the allowed non-renormalizible terms and so avoid suppressing the 
density of coherent S scalars. 

Although we have concentrated on the Z 3 -symmetric NMSSM, the coherently 
oscillating scalar mechanism for eliminating weak scale domain walls in SUSY mod- 
els should have rather general applications to any SUSY model with a weak scale 
domain wall problem, for example, SUSY models with spontaneous CP violation 
or spontaneous R-parity violation at the lOOGeV to ITeV scale. In these cases, 
given the sensitivity of the strong CP parameter 9 to explicit CP violating terms 
and the baryon number violation rate to explicit R-parity breaking terms, it may 
well be advantageous to have a discrete symmetry breaking expectation value which 
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vanishes at zero temperature, as happens in the case of the coherently oscillating 
scalar mechanism. 
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